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Abstract 

The  tensor  virial  equations  of  motion  are  developed.  They 
are  found  to  be  tensor  equations  of  second  rank  which  upon 
contraction  give  the  usual  scalar  virial  equation.  The  tensor 
virial  equations  may  be  apolied  to  anisotropic  systems  in  the  same 
manner  that  the  scalar  virial  equation  is  applied  to  isotropic 
systems.  Several  applications  of  the  tensor  virial  equations 
are  considered.  The  diffusion  of  ions  through  a magnetic  field 
and  the  diffusion  of  molecules  through  a gas  are  calculated. 

The  derivation  of  the  Navier  - Stokes  equations  and  the  Reynolds 
stress  tensor  for  a turbulent  flow  is  developed,  leading  in  a 
natural  wav  to  Prandtl's  mixing  lenrth  ideas.  The  dynamics  of  a 
self  gravitating  oscillating  gas  cloud  Is  investigated.  The 
expressions  for  gravitating  homogeneous  ellipsoidal  regions  are 
worked  out  for  use  in  problems  involving  inverse  square  inter- 
actions, e.g.  galactic  dynamics,  clouds  of  charged  particles  , etc. 
Because  of  the  applicability  to  hydromagnetics,  the  stress  tensor 
formulation  of  the  tensor  virial  equations  is  developed  in  the 
last  section. 
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Introduction 


0.1  The  use  of  the  virial  theorem  in  kinetic  theory  is  well 


known.  It  does  not  seem  to  have  been  remarked  with  sufficient 


generality  that  the  virial  equation  can  be  extended  to  a set  of 


equations  between  tensor  components.  A3  the  formulas  of  the 


first  section  show,  these  relations  are  essentially  equations  of 


motion  for  the  ordinary  moment  of  inertia  tensor  of  an  assembler 


of  particles.  The  usual  scalar  virial  theorem  results  from  the 


tensor  equations  bv  the  process  of  contraction.  The  tensor 


virial  equations  are  useful  in  problems  of  the  dynamics  of  complex 


systems,  especially  in  hydrodynamics  and  transport  theory  where 


the  actual  integration  of  the  equations  of  motion  is  not  practical. 


One  may  then  obtain  approximate  information  by  investigating 


the  behavior  of  the  moment  of  inertia  of  the  system  or  of  a 


subsystem,  e.g.  an  eddy.  The  tensor  form  of  the  equations  permits 


one  to  obtain  results  for  anisotropic  motions  wherever  the  scalar 


virial  can  he  applied  to  isotropic  motions. 


The  method  hes  proved  useful  in  astronhysical  problems. 


particularly  those  dealing  with  the  dynamics  of  a finite  region, 


e.g.  the  dynamics  of  an  individual  star,  clusters  of  stars,  and 


galaxies  '.  We  shall  treat  some  illustrative  examples  in  trans- 


port theory.  We  shall  find  that  the  main  advantage  will  ^e  a 


means  of  sidestepping  .uninterest ing  detail,  such  as  the  precise 


density  distribution,  allowing  one  to  compute  direct1 y the  more 


pertinent  quantities,  e.g.  mean  diffusion  or  expansion  rates, 


mean  pressures  and  viscous  stresses  etc. 


We  shall  derive  the  Navier  - Stokes  equations  from  the  tensor 
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virial  equations  and  show  that  the  parallelism  of  viscous  and 
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^/turbulent  stresses  is  inherent  in  the  theory  without  any 
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ms  a basis  for  the  turbulence  theories 
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to  be  found  elsewhere  ” ..  The  generality  of  our  results  makes  them 


applicable  to  the  construction  of  a statistical  magnetohydro dynamic 

‘ ■*  S ^ s'-  * 

theory- of  the  same  nature  as  existing  turbulence  theories. 


It  is  well  known  that  the  scalar  virial  theorem  holds  in' 


quantum  mechanics*,  an  extension  of  the  tensor  virial  theorem 


to  the  quantum  mechanical  case  should  be  Dossible  , but  will  not 


be  attempted  here. 


"**= 


Lagrangian  Formulation 


1.1  Consider  a system  of  particles  in  a space  with  the  cartesian 

coordinates  x* . We  shall  for  the  moment  fix  our  attention  on 

the  vth  particle.  We  reoresent  its  mass  bv  rri  and  its  nosition 

V 

bv  x1.  We  assume  that  all  velocities  are  small  c imDared  to  the 
v 

speed  of  light.  The  Lagrangian  is  written 


L = T - V 
v v v 


(1) 


where  yT  and  vV  reoresent  the  kinetic  and  potential  energies 


Let  us  introduce  the  svmbol  vQ^  to  denote  those  forces  for 


which  a scalar  ootential  does  not  exist.  In  this  category  we 


include  al]  forces  of  constraint.  The  ^ are  defined  as  the 


differential  coefficients  in  the  expression  for  the  work 


& W = Q.6  . x1 


(2) 


v v i v 

Repeated  italic  indices  imoly  summation  convention:  Greek  do  not. 
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The  &yx  represent  an  arbitrary  displacement  of  the  vth  particle 


under  the  effect  of  all  those  forces  not  included  in  V, 

v 


Lagrange’s  equations  for  the  vth  particle  become 


d /SvL  'N 
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3 x- 
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n 
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MultJnlyinp  by  vx,J  gives  us  the  tensor  equation  of  second  rank 


which  "tfiv  be  written  as 


. i 3 L 
x J v . x 
V “ + v 

V 


3 L 
v 

3 xi 


+ vxVQi 


Consider  now  the  sum  over  a number  of  particles  in  the 
space.  We  shall  indicate  the  summation  by  . In  practice  this 

sum  is  usually  effected  by  summing  over  all  particles  in  a simply 
connected  region  of  space,  though  in  some  diffusion  problems 


(cf 


2.12)  this  is  not  the  case.  We  define  the  tensors 


vo, i 


(6) 


and  the  scalars 


(7) 


(e) 


(9) 


T =21  yT  (10) 

v 

v = 21  V (11) 

v 

Thus,  L,  T,  and  V are  the  Lagrangian,  kinetic  energy,  and  ootential 
Gricr^y  ctively  of  the  particles  over  which  we  are  summing. 

(4)  may  be  rewritten  as 

“dt  Jij=  2Tij  + ^ ij  + -^"ij  (12) 

These  tensors  admit  of  a straight  forward  in ternreta tion . 

The  diagonal  components  of  J, , are  just 
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We  do  not  use  summation  convention  on  Greek  indices,  ^he  diagonal 


terms  of  ^ are,  then,  essentially  the  rate  of  change  of  the 
inertia  tensor.  The  off  diagonal  terms  represent  the  angular 


momentum  of  the  system.  Now  dL/S^x"1'  is  just  ^m^x^.  Thus 


• i 4rl 

vmvx  vx 


The  spur  of  T^  is  simply  T,  the  total  kinetic  energy  of  the 


system.  Hence,  we  call  T.  . the  kinetic  tensor.  QD  , ..  on  the 

ij ^ 


other  hand,  is  given  by 


<£>  - 


.i  3V 


vX  d 


so  that  , if  V is  a homogeneous  function  of  degree  n cf  the  vx  , 


Euler’s  theorem  gives 


^ii  = - nV 


(15) 


i.e.  the  spur  of  is  some  number  --n  times  the  potential 


energy  of  the  system.  Hence  we  call  . the  Potential  tensor 


of  the  system. 


(5)  mav  be  resolved  into  its  sT’Tn,ne  trie  and  ant 3 sTrmmotric 


parts.  We  define 


'iJ  " % 


2 v 


m x1  x^ 

V V V 


K.  . = 1 / m x1  x^  - x ^ 

2 v v v v v J 


h hi 


so  that 
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Further,  we  let 


^--1 

— = § 

l 1 

l 

.4. 

A=§ 
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Thus  we  have  the  symmetric  and  anti-symmetric  equations 

"dtS  Iij  = 2Tij  + Tlj  + Mij 


(20 


(21 


(22 


(23 


(24 


d 

dt 


'ij  ' "ij 


(25 


The  symmetric  equation  gives  the  manner  in  which  the  moment 
of  inertia  tensor,  1^,  varies.  The  antisymmetric  equation  give 
the  rate  of  change  of  the  angular  momentum  of  the  region. 


It  is  well  known  that  if  the  interaction  of  the  ^articles 
mav  be  exoressed  in  terms  of  central  force  fields,  they  do  not 
transfer  angular  momentum.  Hence,  In  the  absence  of  external 
fields  we  have 


■ij  ' Hij  " ° 


(26 


if  Vi*.. 


\ a 


i 


so  that  (14)  and  (15)  may  be  rewritten 


dt2  Iij  2Tij  + +^~i  j 


(27) 


dt  i j 


K.  . = 0 


(28) 


(See1  Appendix  A and  B for  a more  general  development  and  further 


discussion. ) 


1.2  To  obtain  the  conventional  scalar  viriai  equation  we 


contract  (1.2.5)  and  obtain 


~ J = 2T  +<3?  +X1 

dt 


where  in  general  we  define 


X s X.« 


We  note  from  (1.2.8)  that 


J = J.  . 

i i 


d I . . 
“ li 


= 1 d 


2 dt  v 


where 


Using  (2)  we  ma??-  write  in  place  of  (1) 


-^2  I = 2Z  +$  +X1 


Eg1 
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Appli catior s 


2.1  It  will  be  our  ournose  in  this  section  to  illustrate  some  of 
the  methods  of  abdication  of  the  tensor  virial  equations  to 
transport  problems  and  associated  kinetic  nhenomena. 


2.11  Consider  a simply  connected  cloud  of  electrons  and  protons 
in  a uniform  magnetic  field.  We  assume  that  the  net  charge  of 
the  cloud  is  zero  and  recombination  is  negligible.  The  Lagrangian 
for  the  vth  uarticle  is 


vL  = i vVV1  + 4- 


A1  v1 

V V 


where  is  the  vector  potential  reuresenting  the  magnetic  field, 

vq  the  charge.  It  follows  that 

= vravyl  + va1 
o.  v-l  c 


— m ^ rr  A. . A ^ 


9vL  = vM  vVJ  dvA’ 

3VX"  c VJ 


We  obtain  from  (1.1.12),  after  some  rearranging 


ft  ^ v>V\ 


^ 7L*  - 


vmvvlvvJ  - vx^ 

c dt  y ox-1 


If  we  assume  the  magnetic  field  to  be  of  magnitude  B and 
in  the  positive  z direction,  then  we  mav  take 


K = -I  By,  A = 1 Bx,  A = 0,  Mx  = 0 
x 2 7 2 z 3t 

Thms , we  obtain  the  three  equations 


fc  t 


—r  ^ — m v x^  = 
dt  v v xv 


I 1 i v 

•<  m v v-'  + xJ  — 
|V  v xv  v 


q 

— v B 
c v v 


dt  vmvvyvx 


5 j n v - x^  vq  v v B 

^ — <v  v yv  v v x 

V c 


“35-4-  v“vvsvx  = 


vmvvzvv< 


Summing  over  the  whole  cloud  the  off  diagonal  terms  average  out 
because  there  is  no  correlation  between  v^,  and  v^',  or  x*  and  x^. 
(i  ^ j)  as  we  sum  over  v.  We  are  left  with  the  diagonal  terms 


-rr  ^ m v x = 
dt  v v xv 


^fvm(vv. 


)a  + vq  X 


v B~\ 

) V > 


— v v y > 
c J j 


■dt  4-  vmvvyvy  = 


> f m(  v )2  - vq  y v B >• 
— < v v y — " v-^v  x j 


m v z = 
V V zv 


vm(vvz)' 


We  see  from  (9)  that  we  have  free  expansion  in  the  z direction. 


Consider  (7).  We  must  compute 


vqvXvvy  B/c*  We  sha11 


describe  the  average  position  of  a oarticle  between  successive 


collisions  by  means  of  the  center  of  mass  of  its  trajectory.  The 


projection  on  the  xy  plane  of  the  velocity  is 


v =-,/ ( v )2  + ( v )2 
v xy  — v v x v y 


The  radius  of  curvature  of  the  nroiection  on  the  xy  olane  is 

..A  = CWt tt  vM 


We  shall  describe  the  projection  of  the  trajectory  onto  the 
xy  plane  as  shown  in  Figure  1.  The  ft)  axes  are  oriented  so  that 
the  trajectory  is  symmetric  about  the  J axis.  The  rj  component  of 
the  moment  of  the  trajectory  about  § = a is 

1^  (a)  = / ( £ - a)  ds 

where  ds  is  an  element  of  length  along  the  trajectory.  We  have 

S^/2  # 

L (a)  - f ( A co sco  - a)  doo 

71  -Ws 

s ince 

= /A  cos uo  and  ds  = A dca 

where  co  is  measured  from  the  positive  % axis.  We  obtain 

1^  (a)  = A j^2  A s in  ^ - a 

if  ^ = a,  T)  = 0 is  the  center  of  mass  of  the  trajectory,  1^  (a)  = 0 
giving 


a = 


y, 
s In  2 

V 

2 


(12) 


Returning  to  the  xy  coordinate  system,  if  we  take  the  center 
of  mass  of  the  trajectory  to  have  an  abscissa  xQ,  then  the  center 
of  the  trajectory  is  seen  from  Fig.  2 to  be  at  xfl  + a cos<^. 

The  trajectory  is  traversed  in  the  direction  indicated  for 
vq^0,  B^>0.  If  0 is  the  angle  between  the  positive  x axis  and 
the  line  drawn  through  the  center  of  the  trajectory  and  a given 
point  on  the  trajectory,  then  for  the  point  on  the  trajectory 


x = x_  + a cos 


■f  A cos  0 


(13) 


V S=  -V  cos  e 

y xy 


Hence,  using  (12),  (13),  and  (14) 


y. 

: v 


Uf  J 


5! 


(14) 


• -}  / • 


XV  = 

7 


s in 

yAy  O ±11  ^ A p 

-V^^CXqCOS © +/  \ qp^  COS  ^3  COS©  +/\  COfl  0 ) 


(15) 


We  shall  now  average  over  the  tralectorv  operating  with 

~n+<p  + y^/2 

7 

■K+cp  -'i'/2 


We  obtain 


Ixv  )„  — V X COS  CP 

y 0 xy  oT 


sin  ^ a sin2  ^ 2 

TTT  + A -g  003  S* 


(16) 


- A /^  , sin v/,  ^ \ 

2 C1  cos2r y 

where  the  subscript  0 indicates  that  the  average  over  0 has  been 


carried  out. 


We  now  average  over  all  possible  orientations  of  the  projec- 
tion of  the  trajectory  by  operating  with 


h d? 


We  obtain 


iin2^T 


(17) 


Now,  if  the  mean  free  path  of  a particle  is  X,  then  the  probability 
of  its  traveling  a distance  s between  collisions  is  exp  [-s/xj  (l/X). 
But  '•f'  = s//\  . Thus,  the  average  over  s gives 


^XVj^O<jps 


g- ^ l exn(-s/X) 


w 


' v 

r-  -i 


d(f) 


y/\v 


xv 


,x 


2A  ( . / 2A 

— \ du  exp  - — U 


4 2 

sin  u 
u2 


-± 


which  may  be  integrated  to  give 


(xVeVs 


Av 


xx  1 zA 

2 I X 


Arc  tan  In  ML  +^rv 


7^)  J 


-1 


Using  (11)  we  obtain 


^ ..x  ,v_  B = vm  l^jSEL  V 


c v v y 
Thus 


■ X /\  / 'X2  ^ 

Arctan  ^ In  ^H-^g J 


v *— 


»»<„**>  + 


ITvW5]  = 


f 


m . . . . 
v v X v x 


1 (vVx> 


+ I (vVxy)2 


sA  / 


\ 


Arctan  ^ ^ In  (l  + 


However,  we  exuect  the  velocity  to  be  statistically  iso- 
tropic over  the  xy  plane.  Thus 

( y )2  = I ( V )2 
Vx'  2 vv  xy' 


and  we  obtain 


v l_v  v x 


m(  v)^  + ~ ,.x..v_  B 


c v v y 


= f (x/A) 


vm(vvx)‘ 


(18) 


(19) 

-) 

(20J 

-1 

(21) 

(22) 


whe  r e 
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f (x/A)  = 


Arc  tan 


7P/J 


We  note  that 


f (x)  ^ 1 - £ X2  + X4  + 06  (x) 


(x)~  i - 2 + In  (1  + x2)  + + 0~6  (x) 


We  rewrite  (7),  (8),  and  (9)  as 


(83) 


dtX  vmvTxvx  = f(X/A)  .21  vmW 


(24) 


vmvVyv7  = fU/A) 


vra(vV 


vmvVzvZ  =. 


m(  v ) ' 

V V 2 


(24),  (25),  and  (26)  are  not  sufficient  to  determine  the  form 
of  /Q(x,  y,  z),  They  renresent  a restriction  on  the  system  of 
particles  under  consideration,  though  they  by  no  means  determine 
a unique  state  of  the  system.  The^  give  no  information  concerning 
the  form  of  the  spatial  distribution  of  the  particles.  Hence  in 
order  to  use  (24),  (25),  and  (26)  we  must  supply  a form  of  cistri- 
bution  over  space.  (24),  (25),  and  (2G)  will  then  tell  us  how 
the  scale  of  the  distribution  varies  with  time.  With  three 
equations  it  is  obvious  that  the  form  of  the  distribution  may 
be  characterized  b^  three  narameters,  one  parameter  describing 
the  scale  of  the  distribution  in  each  of  the  three  directions. 

At  first  thought  it  seems  a serious  drawback  that  our 
equations  will  not  suorly  us  with  the  form  of  the  spatial  dis- 
tribution as  well  as  its  variation  with  time.  However,  it  must 


14 


be  -emembered  that  the  main  difficult"-  in  oroblems  in  transport 
theory  when  the  solution  is  attempted  frm  the  Boltzmann  equation 
is  the  tremendous  mathematica]  comolexitv  of  the  oroblen.  The 
complexity  is  brought  about  by  our  having  to  determine  simultan- 
eously the  form  of  the  spatial  distribution  as  well  as  the 
variation  with  time.  Hence  we  must  cut  down  on  the  amount  of 
information  that  we  seek  to  obtain  from  our  calculations, 
supplying  the  deficit  by  some  informed  estimates.  The  weakness 
of  our  method  is,  then,  the  source  of  its  main  advantage.  In 
most  cases  the  form  of  the  spatial  distribution  is  not  as  impor- 
tant to  us  as  the  variation  in  time  and  is  usually  known  fairly 
well  anyway.  The  theory  of  stochastic  processes,  for  instance, 
tells  us  that  in  most  cases  the  spatial  distribution  will  be 
gaussian  for  sufficiently  large  values  of  time.  The  end  result 
is  that  in  practice  we  may  assume  a form  for  the  soatial  distri- 
bution characterized  by  three  un5ts  of  scale,  one  in  each 
direction,  and  determine  the  variation  with  time  of  the  three 
units  of  scale  and  ultimately  of  the  spatial  distribution  . 

For  the  problem  at  hand,  any  initial  snatial  distribution 
will  eventually  diffuse  into  a gaussian.  However,  the  puroose 
of  this  oarer  is  to  exhibit  the  general  mehtods  rather  than  to 
solve  accurately  specific  problems.  Therefore,  to  save  computation 
we  shall  assume  that  the  particles  are  constrained  to  a homo- 
geneous distribution  within  a rectangular  parallelepiped  with 
sides  s*y  (t)  and  a,(t)  and  center  at  the  origin.  To  compare  this 
artificial  problem  with  a real  problem  in  nature  we  identify  the 
a (t)  and  a,(t)  as  characteristic  scale;  of  the  distributions, 
say  as  twice  the  mean  deviation  of  the  final  gaussian  distribution. 


1 8 


The  expansion  is  restricted  by  tie  equation  of  continuity 


to 


vx(x,t)  = 

^ xy 

(t) 

X 

(2*7) 

vy(y,t)  =■ 

f]  xy 

(t) 

7 

(28) 

v (z.t)  = 
z 

« 

1 1 z 

(t) 

z 

(28) 

It  follows  that 


d a 

ST  _ 

dt 


h 


xy  (t)  axy 


(30) 


etc. 


Now,  in 


vmyvxvx  all  motions  of  smaller  scale  than  the 


dilatation  average  out,  leaving  us  with 


vTxvx  * JdV  fixy  (t) 

= M ^XV  lt)  ®IT  (t) 


;mvvzvz  = 12  M ^ (t)  az  (t) 


where  M is  the  total  mass  of  the  cloud. 


M = 


m 


(31) 

(32) 

(33) 


' I 


% , 
5*  I 


\ 


We  let  (v  T2  represent  the  mean  square  velocity  in  the  x and  also 

A 


in  the  y direction  (cf.  (21)),  and  (v  ) in  the  z direction.  Then 

z 


vra(vvx)  = MT^T' 


(34) 


etc.  Tv~7  and  (v_5 , it  must  be  remembered,  include  the  small 

A Z 

scale  motions  as  well  as  the  dilatation  velocity.  Hence  the  total 

-2 


energy  density  (v)  /2  is  given  by 


!j 

t I 


7 ; 
I i 
i 
\ 


• -lAjygsL;  


. wnw  c— .-»  - 


■ or‘  - : 


f ’ - r ~V_ . — 
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(35) 


— 2 

Conservation  of  energy  implies  that  v is  independent  of  t. 


Using  (31)  and  (34),  we  may  now  write  (24)  and  (25)  as 


— K 

dt  Ln 


xy(t)  a*y  (t)J  = 12  f (X/A ) TTTC 


(32)  gives  (26)  as 


& u 


z(t>  4 (t)J  = 12  Try2 


Using  (30)  we  obtain 


dt  [axy(t)  £ axy(t)J  = 12  f(X/A)  T rrs 


It  [az(t)  di  az(t) 


)j  =12  T^T2 


(37) 


Let  us  consider  several  special  solutions  of  (36)  and  (37). 
If  B is  sufficiently  small  or  (O  sufficiently  large,  then  X<<  /\ 
and  f (X/A)~  1.  Hence  the  diffusion  is  unhindered  by  B and  we 


have  an  isotropic  adiabatic  expansion.  We  take 


a (t)  = a (t)  = a(t) 

«*■  J 2 


tkt*  = try2  = \ 


1 rrS 


In  place  of  (36)  and  (37)  we  obtain  the  single  equation 


It  [a(t)lt  aCt) ] = 4 vJ 


which  may  be  written 


dSa  . /da  \ 

^ (dt  j 


. —2 
= 4 v 


(38) 


\ 
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It  should  be  noted  that  sufficient  expansion  will  eventually 
increase  X to  where  it  is  no  longer  less  than  /\  . 

It  we  consider  the  other  extreme  of  a stronp  field  and/or  a 
low  density,  X^A,  and  The  cloud  expands  only 

in  the  z direction.  The  enerry  of  the  expansion  is 

1*  M 

z 


i Xdv/°  [hz<t)  z]2  ft) 


(39) 


The  energy  of  all  notions  of  smaller  scale  than  the  parallelepiped 
is,  then, 

, „ / da_  2 

(40) 


I"  ^ 


M_  fe) 
" 24  \ dt  / 


(We  note  the  division  of  the  motion  into  two  components  defending 
upon  the  scale  relative  to  the  region  considered.  More  will  be 
said  on  this  in  2.21J  Assuming  this  energy  to  be  distributed 


isotropically  over  the  three  dimensions,  we  have  T as  just  the 

z z 


energy  in  the  dilatation  velocity  olus  one  thir 
in  the  small  scale  motion. 


d of  the  energy 


/da  \2  , 

. JL 

/da  \2 

f __z  + I 

- M v2 

—5.  ) 

\ dt  / 3 

_2 

24 

^ dt  ) 

M 

6 


T2  + -'  ( ■ 


■ feT 


(41) 


Thus , (o7j  becomes 


,2q 

u a 


* dt2 


da. 


dt 


= -4  v2 


(42) 


New  let  us  put  no  restrictions  on  X but  constrain  the  particles 
so  that  there  is  no  expansion  in  the  z direction.  This  can  be 
done  by  considering  the  region  confined  between  two  material 
planes  < rpendicular  to  the  z axis. 


f i 

3 * 


it 


: fH ir'WMHIttiHI  ’ *ra 


'•  y^iem 
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Analogous  to  (39)  , it  can  be  shovn  that  the  energy  of  the 
expansion  in  the  x and  y directions  is  each  given  by  (M/24 ) (da^^/dt ) 
Thus,  the  motions  of  smaller  scale  than  the  region  have  an  energy 

M /da 

1 M ~2  _ JL  —Z3.  ) (43) 

2 v 12  V dt  / K 1 

And,  assuming  isotropy,  T and  T are  just  equal  to  the  energy 

xx  yy 

of  the  dilatation  plus  one  third  of  the  small  scale  energy  given 

in  (43).  - -] 

T = T = JL  &SZ  f + M i j2  .i  M 

Xxx  yy  24  \ dt  / 3 2 12  \ dt  / 

(44) 


■ h~'  * * ft?;. 


(36)  becomes 


fi-im/A)]  (%)  = 4f(vA).vS 

(38),  (42),  and  (45)  are  all  of  the  same  general  form 


(45) 


y ix  + 
y . 2 
dx 


a (&f  = 

\dx/ 


where  a and  P are  constants  if  we  neglect  the  variation  of  X with 
the  expansion.  To  integrate  (46)  we  let  p = dy/dx  and  write 


±JL  = d£  _ d£ 
,2  dx  p dy 


(47) 


W8  obtain 


(§SLf  = £ 
\ dx  / a 


, C 
1 " 7** 


(48) 


where  C is  the  c istant  of  integration. 

Let  us  assume  that  when  t = 0 there  is  no  expansion  and 
y = yQ.  (48)  becomes 


[49 


I 


**»• 


„ hr, 

m 


I: 

V ■' 


1 


N 


= 

dx 


-ft) 


Thus,  we  obtain  the  solution  by  quadrature, 


2a  2a 

y - y 

c.  y 


For  (38)  we  have  a = 1 and  6 = 4 v^.  (49)  and  (50)  give  the 


rate  of  expansion  of  each  face  of  the  parallelepiped  as 


1 d 


2 dt 


.(t)  = V v' 


a(o) 

[_  a ( t ) _ 


2 1/2 


(51) 


and  the  length  of  each  side  as 


a2(t)  - a2(o)  = 4 v2  t2 


For  large  values  of  t we  note  that  the  velocity  of  expansion 


1 d 


2 dt 


a(  t)‘ 


(53) 


in  agreement  with  the  usual  result  for  expansion  into  a vacuum. 


For  (42),  a = 1/3,  0=4  v^.  (49)  and  (50)  give 


az(c)  ]2/Z{  ^ 

^TtT  f 


t - [a2/3(t)  - a2/3(o)  ] 1/2  [j^t>  + 8J>V>]  (6 


2 V 3 v2 


For  large  values  of  t,  the  rate  of  expansion  is 


t dar 


2 dt 


■/3  "v2 


Finally,  for  (46) 


a = 1 - — f(X,‘A) 


(3  - 4f(X,/\)  v2 


itvar.^ 


T-—  ■ 


K\ 


(56) 


**  r Wib  <%  ■SbbI  (4- 
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We  obtain  the  diffusion  rate  as 


1 daxy  = f(X/A.)  v2 

2 dt  1 ” 3 f (>/A) 

For  large  values  of  t. 


1-/2 


I \ 1 2( 

a _( t) 
xy  J 


(l-f  (X/A)/3  / 1/2 


1 daxv 

2 dt 


_ ~ 1/2 

f(X/A)  v2 

1 - 4 f (x/A  ) 


We  note  that  if  X<</\,  f(X/A)~l  and 


(57) 


(58) 


, da 

1 ...iHE 

2 dt 


3 v^ 
2 


(50)  results  in  an  integral  which  must  be  evaluated  numerically 
except  for  special  values  of  f(x/A). 

If  one  wishes  to  consider  further  complications  such  as  a 
uniform  distribution  of  neutral  particles  which  inhibit  the  dif- 
fusion of  the  charged  particles,  one  introduces  an  additional 


force  into 


F^-  x^  due  to  collisions  of  the  charged  with  the 


neutral  particles, 


2.12  As  an  examnle  of  a diffusion  problem  consider  a snace 

filled  with  a homogeneous  distribution  of  neutral  particles.  At 

time  t we  mark  every  particle  within  the  rectangle  with  sides 

a « a , and  a and  center  at  the  origin.  We  ask  how  these  marked 
x’  y ’ z 

particles  will  be  spread  out  at  some  subsequent  time  t.  As  has 
been  pointed  out,  the  tensor  virial  equations  do  not  determine 
the  form  of  the  spatial  distribution,  but,  upon  assuming  some 
spatial  distribution,  indicate  how  it  will  vary  with  time.  It  is 
convenient  to  express  this  idea  by  stating  that  we  use  the 
tensor  virial  equations  to  determine  the  dynamical  conditions  for 


a given  scale  without  inquiring  into  fluctuations  of  smaller 

scale.  Now,  it  can  be  shown  from  the  theory  of  stochastic 

processes  that  the  asymptotic  distribution  as  t — *od  is  a gauss ian 

centex ing  about  the  origin.  Obviously  our  initial  step  function 

distribution  rounds  its  corners  and  spreads  out  with  increasing 

time.  It  is  not  our  purpose  here  to  go  into  so  much  detail. 

Thus,  we  describe  the  distribution  by  the  three  lengths  a (t), 

a (t) , and  a (t).  For  t = t we  identify  them  with  a , a , and 
y • * z o x7  y 7 

a . They  correspond  to  approximately  twice  the  standard  deviation 
z 

of  the  distribution. 

The  tensor  virial  equations  reduce  ", o 


_d_ 

dt 


ai(t) 


_d_ 

dt 


a.i 


0 


(1) 


since  the  pressure  on  the  boundaries  just  balances  the  kinetic 
tensor.  The  solution  of  (1)  is  readily  shown  to  be 


ai  ( -t ) = ai(tQ) 


(2) 


_d_ 

dt 


ai(  t) 


(3) 


■■u 


!? 


We  must  now  evaluate  t • The  rate  of  increase  of  a0(t)  at 

o P 

each  face  of  the  parallelepiped  is  (1/2) (da^/dt) . At  time  tQ 
the  expansion  occurs  only  within  a distance  L of  each  face, 
where  L is  the  mean  free  path;  the  particles  farther  into  the 
parallelepiped  are  as  vet  undisturbed.  Cf  all  the  particles  in 
the  parallelepiped,  a fraction  L/  [ap(to)/2]  are  contained  in 
the  slab  of  thickness  L ana  normal  to  the  p direction.  Half  of 


atmpyi  iar.u»* 
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these  particles  are  moving  outward  across  the  lace  with  a velocity 

| . 

Vta)*  • Therefore,  the  fraction  L/a.p(to)  of  the  particles  have 

I 

a velocity  -Jfy*)*  outward  across  a face  normal  to  the  6 direction. 

We  write  that  the  rate  of  change  of  the  characteristic  length 

aft(t)/2  is 

r' 

„ r<r  (t)i  T / 

« rV-J  - ^ (4) 

1 

of  time  t . Putting  t = t in  (3)  and  comnaring  with  (4).  we 
o ■*  o 7 

I 

obtain 

2 

a« 

t^  = % r-  (5) 

4L  -j  (vfi)  2 

1 

Thus,  from  (2)  and  (3)  we  obtain 

ai(t)  = 2 |l  -/  (vp)-<i  tj  ^ (6) 

T 2 /£ 

I 

j Ly  y_ 

dt  ai^  “ |_  t j ^7- 

One  advantage  of  this  method  is  obviously  the  doing  away 

m 

with  the  infinite  diffusion  rate  obtained  at  t = t if  the  oroblem 

o 

I 

were  fomulated  in  the  conventional  manner  in  terms  of  . The 

I 

thermal  diffusion  equation  is  altered  to  give  a finite  rate  of 

■ 

y 

A ■ 

propagation  of  thermal  disturbances. 

1 

2.2  Let  us  use  (1.1*12)  to  derive  the  equations  of  motion  of 

9 

a finite  region  of  fluid.  We  consider  a finite  region  moving 

"W 

with  the  fluid.  We  denote  the  center  of  mass  of  the  region  by 

I 

the  cartesian  coordinates  x^.  We  shall  locate  points  within  the 

r 

«' 

region  relative  to  the  center  of  mass  of  the  region  by  the 

coordinates  §*,  so  that  the  coordinates  of  the  vth  particle 
become 

become 


TSflMIP 


I 

i 

I 

£ 


i 

f 

l 


i = -i  ^ t i 

V V V- 


x'  = -.*■  ^ VJ 

The  velocity  of  the  vth  particle  is 


23 

(1) 


i _ V. 

v v at 


(2) 


which  we  shall  find  convenient  to  write  as 


vvl  = ^ + vui 


(3) 


Cle  arly 


SI  vVx  = 

V V 


v”vU  = 0 


(4) 


We  refer  to  the  vu^  as  the  lo  cal  and  v*  as  the  translocal 


velocity  field  of  the  region.  Thus,  the  local  field  is  the  portion 
of  the  velocity  composed  of  fluctuations  of  smaller  scale  than 
the  region.  Or,  in  other  words,  the  local  field  is  the  portion 
which  does  average  out/  the  translocal  field  is  the  oortion  of 
the  velocity  field  which  does  not  average  out  over  the  region. 

Prom  (1.1.5)  - (1.1.7)  and  (4)  it  follows  that 


Jij*  = MxV  + 


m u^ 
v v v 


V 


i pj  = 


'j  + 


pj 

V V 


2T^j  = Mv^vj  + .^m.^u^^uj 


(1.1.12).  becomes  after  some  rearranging  of  terms 


u1  u j + / j 

V V V>  V 


(5) 


_d_ 

dt 


J 


m £ J u 
v Vs  v 


j 
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i - 


A 


^4 


•'4 


m 


But,  from  Newton’s  equations  we  have 


M — * = 

“ At  +■ 


Thus 


vV  v 


f1  J = 

* \l 


vmv 

V 


u*  u^  + 


jl  pj 

V 


for  the  region,  independently  of  the  motion  of  its  center  of  mass. 
We  shall  now  use  (7)  to  evaluate  .fK  We  decompose 


F into  two  portions. 


pi  = pi 


vF^  represents  the  mean  external  force  e.g.  gravitational  or 


electric  fields.  We  shall  call  it  the  translocal  force  field. 


Finally  yf  denotes  the  short  range  forces  of  collision  of  the 


molecules . 


Since  F ^ is  indenendent  of  oosition  in  the  region,  it 
vo  ’ 


being  the  average  over  the  region  and  at  most  a function  of  the 


tyne  of  particle,  we  have 


F 3 f1  = 0 

v o v 


from  (4).  (7)  becomes 


vravu‘Jv 


A .J ) * 


Vs- v v v I 


A ri 


5 f' 
V3  V 


(10) 


Now,  the  f , being  collisional  forces,  cancel  out  over  the 


interior  of  the  region.  The  only  contribution  can  come  from  the 


surface  of  the  region  where  collisions  take  place  with  particles 


net  in  the  region,  so  that  the  equal  and  opposite  force  on  the 


other  colliding  particle  is  not  included  in 


. Thus,  we  find 


it  convenient  to  define  the  usual  stress  tensor such  that 


m 
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l 

i 

i 


» 

r 

k 


! 

* 

i 

> 


I 


| 

» 


t 


\ 


i 


i C 

d dSa  is  the  force  in  the  i direction  across  an  element  of  area 
dSa  normal  to  the  a direction,  6*°’  represents  the  force  exerted 
by  the  matter  on  the  nositive  side  of  dSa  on  the  matter  on  the 

5 

negative  side.  This  is  the  customary  definition  in  elasticity 

and  electromagnetic  theory.  We  shall  reolace  by  integration 

v 

over  the  elements  of  volume  dV. 

The  introduction  of  a stress  tensor  and  the  associated 
processes  and  parameters  such  as  integration,  differentiation, 
pressure,  viscosity,  etc.  requires  a limited  form  of  continuity 
in  our  hitherto  unrestricted  system.  Our  notion  of  infinitesimal 
becomes  that  of  the  physical  infinitesimal,  viz.  that  the  smallest 
elements  of  volume  dV  that  we  consider  must  be  sufficiently 
large  as  to  contain  a large  number  of  particles.  If  n reoresents 
the  number  of  particles  in  dV,  then  the  fluctuation  of  n is  of 
the  order  of  ~V^rT  . We  must  require  that 

-v/n~  < < n 

in  order  that  our  averaging  process  over  dV  have  a smoothing 
effect.  To  be  treated  as  an  infinitesimal,  dV,  of  course,  must 
be  of  smaller  scale  than  the  ohenomena  in  which  we  are  interested.  | 

The  alternative  to  these  restrictions  on  dV  is  to  consider  an 
ensemble  of  systems  so  that  the  average  may  be  carried  out  over 
the  given  dV  in  each  member  system  rather  than  .iust  a single 
system. 

We  rewrite  (10)  as 

j dV^>  s V)  = JdV(^uM)  + j dSk5i<^k  (11) 

Gauss'  theorem  gives 


I 

! 
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{^|dV*j!l  = j dV u1  u “ ) + j dV-^k(| (12) 


= jdV(/0uV)  +^n 

3| 


(13) 


-3djk/a|k  reoresents  the  force  in  the  1 direction  oer  unit 
volume.  The  translocal  portion  of  - a^Vaf1*  given  to  acceler- 
ating the  entire  region,  cannot  contribute  to  fav^dd^/df  k. 
This  is  readily  shown  if  we  note  that  the  translocal  portion  of 
-3d^k/^ik  is  defined  as 


agJk 

^s=k 


= ^>€3 

at 


Then 


_ * dvJ’  ( 

" dt  J 


&v/0 


by  (4). 

The  remaining  portion  of  3c5^k/3fk  is  local  to  the  region 
and  is  responsible  for  changing  the  shaue  of  the  region.  Since 
6 ^ is  the  only  stress  field  present,  it  follows  that  if  the 
local  component  of  3d>  /3§  vanishes,  then  the  region  is  in 
equilibrium.  We  have  dJ^/dt  = 0 or. 


and  (13)  becomes 


it  $ dV/°s 3ui  = 0 


= \ dV^,  uiu<1  + \ dVd1^ 
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.3 


which  has  the  solution' 


= - /O  u u 


(16) 


Prom  (8)  and  (l7)  we  obtain 


= -2TJ 


F -■  + ( dS , <5ik 

> o J k 


(6)  becomes 


dt 


L_  P 1 - ( 

7 v 0 J 


dV-~(/oUluk) 


■j 


,TT  3 i i k-v 
dV- — ^ (iO  u u ) 

df 


(13) 


Consider  now  (^u'hi^).  If  there  is  no  shearing  of  the 

region,  we  expect  no  correlation  over  the  region  between  and 

* • • 

u^.  Hence  (<juu^)  Is  zero  if  i ^ j.  The  diagonal  terms  are 
nonzero  and  give  the  pressure.  If,  however,  there  is  an  overall 
shearing  of  the  region,  we  see  that  -there  will  be  a correlation 
between  u*  and  u^  for  i j.  Let  the  shearing  be  represented  by 


"We  see,  then,  that  ST1''  represents  the  integral  over 
the  region  of- the  Reynolds  or  kinetic  .stresses  , wbi  be 

represented  bv  pressure,  2TCTa  , and  shearing  stress  -2TaP. 

, as  In  (1.1.12),  represent  the  integral  over  the  region 
of  all  other  stress  fields,  e.g.  gravity,  electric  fields,  etc. 
If  we  denote  by  Y^-1  the  sum  of  the  Reynolds  and  all  other  stress 
fields,  then  (1.1.12)  may  be  written 

•yr  J.  . = Y.  . 
dt  i j i j 


If  Yjj  vanishes,  so  also  aoes  dJ-r-t/dt, 
at  greater  length  in  3.1. 


This  will  be  discussed 


’ «**«**»*  • 
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3v  /Sx**  If  the  average  scale  in  the  P direction  of  the  local 

i ^ ft 

velocity  field  u is  denoted  by  Lr  , then  with  a velocity  u'  we 

expect  to  find  associated  a velocity  -L^3v°'/3x^  in  the  a direction. 


Hence,  we  write 


/ a Sv  / fiT8'  3va  / aTa.\  3y 

(pvi  uH)  = -(/ou  LM ) •—  - ( /O  u L ) 


(19) 


for  a p.  We  define 


i / aTcu 

M-  = Uu  L ) 


and  write  more  generally 


(/ouau^)  = -^aa6 


aacP  PP  3v  aa  3y^ 

1 °a  “ ^ .-6  ' ^ _“0 

3xK  3x 


* 


Thus,  we  rewrite  (18)  as 


Fa  +_  ( dV*~v  (5 aa  + 

V O 1 -f—CL 


r 3 PP. 

p 3x1  axp 


*h  f-M) 


which  is  the  equation  of  motion  for  a finite  region  of  fluid.' 


\y\ 


2.21  Let  us  consider  some  snecial  cases  of  the  equations  of 
motion  of  a finite  volume  of  fluid  expressed  in  (2.2.22).  If 
we  consider  an  incompressible  flow,  then 

= o 

3X1 


and  the  last  term  in  (2.2.22)  becomes 


3 /»%  351)  = ( dV  Si1  ± 
dx1  V 3x/  ) axa  ax; 


If  there  are  no  compressibility  effects,  then  to  a fair  degree  of 


29 


approximation  the  variation  of  ji  in  the  direction  of  the  flow 
may  be  neglected.  Hence,  we  are  able  to  rewrite  (2.2.22)  as 


a 


F a + \dV  — aa 
v o 


fdV 

J (5  3x,J 


BP  3s“*\ 

axP  J 


We  note  that 


= - P 


where  p is  the  pressure  in  the  cl  direction.  Thus,  if  there  is 
an  approach  to  isotropy  in  the  local  motions,  we  write 

p = - - 6\ i (3: 


= y-(/oUiu.) 


and  the  second  term  in  (2)  becomes 


dV  ^Za.6aa  = - dV  ^ 

dSE*  3xa 


It  Is  of  interest  to  compare  (4)  with  Batchelor’s  valueb  computed 
by  B'ourier  transform  methods  from  the  theory  of  isotropic  turbulence, 
Again,  if  there  is  an  approach  to  isotropy  in  the  local 
motions,  we  define 

u = (V?^T)  (6) 


L = 


(V  LiL1) 


so  that  we  may  write 


-7==T  U 

VT 


-J~tT 


**“**■'  • 


30 


Assuming  no  correlation  between  ua  ana  La,  we  obtain 

aa.,  1 ~ T 

= -/o  uL 


(10) 


A • 


P." 


~ p 


(11) 


and  rewrite  (2)  for  the  case  of  isotropic  local  motions 


as 


M—  s 
dt 


a 


v o 


3x 


dx~ 


(12) 


If  we  let  the  size  of  the  region  approach  zero,  we  obtain 
the  usual  Navier  - Stokes  equations  for  an  incompressible  fluid, 


(13) 


dva  _ pa  1.  dp  , 8 

dt  " > axa  ax1 

where  Fa  is  the  force  ner  unit  mass  and  v is  \i/^D  . p and  p 
are  due  only  to  the  thermal  motions,  i.e.  motions  local  to  a 
scale  of  zero. 


2.22  It  should  be  noted  that  aap  and  p*10  or  i 6aa  I , and  therefore 
p and  p are  monotonical  y Increasing  functions  of  the  size  of 
the  r egion  considered.  v°'  and  Its  space  and  time  derivatives,  on 
the  other  hand,  are  monotonically  decreasing  functions  of  the 
size  of  the  region. 

It  is  of  interest  to  note  how  with  the  tensor  virial 

equation  the  velocity  field  v1  naturally  resolves  itself  into 

the  translocal  and  local  velocity  fields  v^  and  u^.  We  note  that 

only  the  statistical  characteristics  of  the  local  field  appear  in 

the  field  equations  of  the  translocal  field,  (2.2.22)  Thus,  we 

are  not  surprised  to  see  the  emergence  of  Prandtl's  mixing  length 

i 1 

concept  in  computing  ( ^ou  U")  for  i / j.  Altogether,  then,  it 


I 


• ii*" 

~S 

-S 


' ! 


i ' 
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would  seem  that  the  statistical  representation  of  Prandtl's 

mixing  length  theory  and  the  more  extensive  treatment  carried 

7 R 9 

out  in  Heisenberg’s  field  equation  * ' and  elsewhere  follow  quite 
naturally  from  the  tensor  viriai  equations  of  motion. 


£9 

] 

3 


2.3  The  tensor  viriai  equations  are  particularly  suited,  to  v 

investigation  of  the  dynamics  of  finite  regions  of  matter.  The  ( 

following  examples  are  typical  of  those  encountered  in  astro-  '? 

5 

physics  in  the  treatment  of  gas  clouds,  star  clusters,  etc.  and 
in  diffusion  problems.  In  the  treatment  of  the  dynamics  of  ras 
clouds  the  tensor  viriai  equations  are  convenient  because  even 
with  very  artificial  constraints  the  question  of  whether  a term 
is  to  be  included  as  a potential  energy  or  as  a kinetic  energy 
is  easily  decided. 


2.31  Consider  the  radial  oscillation  of  a gas  cloud  in  its  own 
gravitational  field.  The  tensor  viriai  equations  do  not  give  us 
the  radial  distribution  matter.  To  prevent  the  mathematics  from 
becoming  too  complex  we  introduce  the  constraint  that  the  density 
be  homogeneous  within  the  sphere  of  radius  r^Ct)  about  the 
origin  and  zero  outside.  Thus,  the  radial  velocity  at  any  point 
within  the  cloud  is 


(r)  = 5 


(1) 


The  motion  of  the  system  is  thus 

parameter  r . Hence,  we  use  the 
o 

It  is  readily  shown  that 

V 


describable  by  the  single 
scalar  viriai  equation  (1.2.3). 


(2) 


i • 

i 


r 


( 


j 
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I = ^Hr./  <3) 

where  M is  the  mass  of  the  sphere.  The  contribution  of  the  radial 


oscillations  to  T is 


T 

L1  10  n\  dt  J 


Let  us  assume  that  the  motions  of  smaller  scale  than  r follow  a 

o 


polytrope  law  so  that 


where 


. TTys/^oj4-1 


r = r ( t ) . 
oo  o o 


A corresponds  to  the  usual  y~  written  when  only  thermal  motions 
are  considered.  For  instance,  if  there  is  no  dissipation, 

A = 5/3  for  turbulence  and  for  the  thermal  motions  in  monatomic 


gases.  If  more  than  one  type  of  local  motion  is  present,  we  write 

/r  \3  (A  - 1) 

2Tg  = m nrr2  ) a (5) 


(1.2.3)  becomes 


m2  r 3(Aa-1) 
a o 00 

3A  -2 

ro  £ 


after  cancelling  out  a factor  of  3M/lu.  Upon  integration  we 


obtain  the  energy  equation 


- 

" 10\ dt/ 


^Jo  /ro< 

3‘Aa  - lTl  ro 


r fV11 

OO  1 


where  c is  the  total  energy  per  unit  mas3.  We  obtain  r (t)  by 
quadrature  as 


M if  ^ f f i‘  »**  * j.l  «-  '*■»(  ■ : '’  -j-i'  ■&-  '■  - ' t»4-1  a <*•>  *1* ^4 C.**  ■ '■&*  i >A  : H.’i  , • £ 9- 


t - t 


= /r° 


dr. 


. r 
j oo 


- 1) 


2GM 


1/2 


33 

(8) 
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2.32  One  finds  in  his  applications  of  the  tensor  virial  equations 

that  he  usually  uses  rectangular  parallelepioedal  or  ellipsoidal 

regions.  In  the  case  that  the  region  is  an  isolated  one, the 

parallelepiped  becomes  inconvenient  because  of  the  extreme 

complexity  of  the  force  field  exoressions  from  which  the  potential 

12 

tensor  is  computed.  The  expressions  for  the  contribution  to 
the  potential  tensor  of  the  mutual  inverse  square  interactions  of 
the  particles  composing  a.  homogeneous  ellipsoid  are  recorded  here. 
If  external  forces  are  Dresent  there  will,  of  course,  be  an 
additional  contribution  to  the  potential  tensor.  The  exoressions 
derived- In  this  section  have  been  applied  elsewhere  to  problems 
In  galactic  dynamics' 


2 


Consider  particles  of  matter  distributed  uniformly  throughout 

a 

the  interior  of  an  ellipsoid  with  semi-axes  a.  Without  loss  of 
generality  we  orient  our  axes  so  that 

(1) 


1>  2>  3 

a ^ a a 


12 


It  can  be  shown  for  gravitational  Interaction  that 

.aa 


= - ~ GM2(aa)2  Na 


(2) 


i 


where  M is  the  mass  of  the  ellipsoid  and 

2 


N 


[(ab^'-'u3)8]'^  j k2  ' E<“’k)] 


N2  = 


[(a1)2  - <a3)2r/a 


E(cj  .k) 1 

k2(l  -k2)  (1-  k2) 


3TV  CMV 

dnv 


(3) 


y u> 

k 


F(uj,k)  and  E(co,k)  are  Legendre’s  elliptic  integrals  of  the  first 

and  second  kinds  respectively. 

— u. — s? 

If  (v  ) represents  the  mean  velocity  in  the  a direction, 
then  for  equilibrium  of  the  ellipsoid,  (1.1.12)  gives 

/ aN  2 3 , cu2  wa 

(v  ) = T-r  GM  (a  ) N (8) 


If  the  ras  cloud  has  rotational  symmetry  e.g.  a spheroidal 
galaxy,  and 


1 ' 
a = 

2 

a > 

3 

a , 

(9) 

then  the  <3>aa 

reduce 

to 

&11  = 

-p,2  2 

= = - 

3GM2 

5 a 

^‘YS) 

(10) 

<£>35  = - 

3GM2 

5a 

(11) 

where 

r- 

2(1 

. y2)3/2l 

s in-^ 

/ 1 2x]t£  m 2x]/2 

(1-y  )J/  - y (1  -y  J 

(32) 

,„3,2  3GM  -33  (g) 

" ' 5a1  ^ Wl 


(13) 


f CL  2 

wouj-u  ue  wo  xx  uo  rcemohasize  the  fact  that  the  (v  ) 


7*4-  l i 4- 

XO  VVV-/L4.J.W  U “ 


are  the  mean  square  values  of  the  local  velocity  field  and  need 
not  he  random  in  the  usual  sense.  We  onlv  require  that  the 
average  >j£  \a  over  the  region  vanish.  Thus  the  rotational  velocity 
of  the  region  as  well  as  smaller  scale  turbulence  and  thermaj 
velocities  are  oart  of  the  local  velocity  field  and  make  uo 


CL  o 

(v  ) . For  ihstance,  if  we  consider  a sohericdal  ralaxy  in  which 
the  velocities  local  to  the  rotation  are  isotronic,  it  fellows 
that  the  rotation  velocity  is 

2 


(O  = 2 


_ (v1)2  - (v3)2 


6GM  vTXl/a3 
— T SK  ("T 
5a1  L U > 


] 

i-  ^(4 

\ a /j 


from  (14)  and  (15). 


Stress  Tensor  Formulation 

3.1  Having  discussed  the  tensor  virial  equations  in  terms  of 

the  kinetic  and  potential  tensors,  let  us  now  state  the  relations 

in  terms  of  stress  tensors.  This  formulation  is  Particularly 

13 

useful  when  dealing  with  problems  in  magne tohydrodynamics . We 
shall  assume  ponderamotive  forces  representable  bv  the  general 
stress  tensor  Then  (1.1.12)  mav  be  written'*1* 

Oj-  ^dV^vM  = J dV/oviv^  + ^ dSkxj6ik  + j dVx^  Zik  (l) 

in  cartesian  coordinates.  It  should  be  noted  that  we  cannot 


See  formulation  from  Newton's  equations  in  Aooendix  A 


36 


express  the  contribution  of  Z1^  to  th  ^otentr'  . ensor  as  • 
surface  integral.  The  net  force  or  the  region  go  r be  expressed 
as  y'dS^Z ^ , but  in  computing  the  potential  tenso^  the  position  at 
which  the  force  is  exerted  on  the  matter  field  is  important 
because  of  J * under  the  integral.  Hence  we  have  action  at  a 
distance  appearing  explicitly.  Using  (2.2.1),  (2.2„3) » and 
(2.2.4)  we  may  rewrite  (l)  as 

M ff1'  = /dsylk  + /av  -4%  zik 


+ 


- /dV^ou1^  + /dV^  uV  + /dS^te3^ 


+ /dV|J  zlk  j 

where  the  comma  indicates  differentiation.  From  Newton's 
equations 

M^1  = /dSk«5lk  + zik) 

Thus,  we  have  for  the  local  motion  the  tensor  equation  valid  in 
any  coordinate  system. 


^ /dV*  u^q*)  = /dV^uV  + /dSkf^(q n)dlk 
+ /dV$Hqn)  Z^ 


(2) 


from  which  we  see  that  the  tensor  virial  equations  for  the 
exterior  velocity  of  a given  region  are  independent  of  the 
interior  velocity  field  and  its  time  derivatives . This  was  noted 
in  passing  in  a more  restricted  case  in  2.2. 

We  note  that 


37 


/dSkfkik  = /d7<? 1 ' + /d"|Jdjk 
Using  (3)  and  (2.2.16),  (2)  nav  be  rewritten  as 


/dV^u^U11)  = /dV|j(qn) 


~ /dV^  u1!  J (qn ) = /dSkf  3(qn)  (6ik  + Zik) 


- /dV (6i^  + 

If  we  are  considering  an  isolated  system,  the  surface  integral 
vanishes  and 

^ /dV^u1! J(q**)  = - /dV^-Uz1^  (( 

—2.^  is  ju3t  the  energy  per  unit  volume  of  the  stress  field 

* * 

Z1J  . Thus,  if  we  contract  (6),  we  have  from  1.2  and  (2.2.17) 
that 

4r  = 2T  + S (I 

at 

where  T is  the  total  kinetic  energy  of  the  local  motions  and 

S = -Z?  (e 


The  symmetric  and  anti-symmetric  parts  of  (5)  may  be  'written 
as  in  1.1  to  give 

% 

~ell}  = 2 /dSk  [? J<qS)  Wik+zlk)  - 1(<j: n)»Jk  + zJ*j] 


- /dV(<51J  + Z1^) 


~ =|  /dSk  |J(qn)(ciik  + Zlk)  - f1(qn)(^k  + ZJk)  (10) 
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where 


IiJ  = A /dV/ofi(qn)  fHqn) 


(11) 


Kij  = | /dV^  |i(qn)  - f*(qn)  U1  (12) 

For  an  isolated  system  the  surface  integrals  vanish  giving 

I1 3 = - /dV  (6ij*  + Z^)  (13) 

dt^ 


^rK1 J = 0 
dt 

Contracting  these  two  equations,  we  obtain 


d2I 

dt 


= 2T  + S 


(14) 


(15) 


It  KS0 


(16) 


3.11  In  the  previous  section  we  have  set  uo  the  tensor  virial 
equations  in  terms  of  a general  stress  tensor  Z*^.  It  would  be 
well  if  we  note  briefly  the  form  of  Z^.for  various  fields.  It 

9 

is  readily  shown  that  for  electromagnetic  fields  the  comnonents 
of  the  three  dimensional  Maxwell  stress  tensor  are  given  by 


= irB**  + d 


^ -g^S 


(1) 


S " t Sij  ^n_B  + D (<0 

in  MKSQ  units.  g^  is  the  metric  tensor.  S is  the  energy  density 
of  the  field. 

In  most  nroblems  in  magnetohydrodynamics  it  is  oossible  to 
entirely  neglect  the  disolacement  current^.  The  electromagnet 1 c 


>-  iHUW 


field  becomes,  then,  a magnetic  field  and  is  a ^i..-e  stress  field 
with  no  associated  inertia  to  be  included  in  oi;r  eqxjations.  ’»ve 
droc  the  terms  involving  the  electric  field  so  tnat 

Z1 j = - g^S  (3] 


S = | glJ  HiB*^  (4) 

Further  simplification  may  be  carried  out  since  the  oermeahility 
of  the  conducting  fluid  usually  aoproximates  very  closely  to 
empty  space  and  is  in  any  ease  isotropic  and  homogeneous. 

In  the  case  of  a gravitational  field  we  have  the  field 
intensity  given  in  terms  of  a scalar  according  to 


so  that 


g1J  b',i,3  = 4’'G/- 


^4 


The  force  oer  unit  volume  Is 


Thus,  In  order  th&t 


= iw  III1  U -J  «J  . 

^ . 1 


F1  = Z\k 


the  s tress  tensor  is  defined  as 


— 

2i;J  _ L_  gljw_ 


do) 
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where 


w = | O'  1 <L> 


(n) 


In  order  to  derive  (8)  from  (9)  and  (10)  we  note  from  (4)  that 
the  curl  of  is  zero. 

The  total  energy  of  the  field  is 

S = -Z* 


W 

4nG 


87tG  ^ ^i 


(12) 


Conclusion 

4.1  It  was  Dointed  out  in  2.22  that  the  ideas  basic  to  Prandtl's 
and  Heisenberg's  formulations  of  turbulence  theory  arise  as  a 
natural  consequence  of  the  tensor  virial  formulation  of  the 
equations  of  motion.  Let  us  summarize  the  more  important  of  these 
principles . 

The  effect  of  multiplying  by  before  summing  over  the 
particles  in  the  region  is  to  cancel  out  the  effects  of  force  and 
velocity  fluctuations  of  larger  scale  than  the  region  considered. 

The  result  is  a separation  of  all  effects  into  the  loco!  and  trans- 
local components.  It  was  then  shown  in  3.1  that  the  local  tensor 
virial  equations  are  independent  of  the  translocsl  field.  (2.2.10) 
shows  that  the  translocsl  field  depends  only  on  the  statistical 
properties  of  the  local  field.  We  found  in  2.2  that  the  statistical 


properties  of  the  local  field  relevant  to  a calculation  of  the 


translocal  field  nay  be  constructed  in  a natural  wav  from  a 
characteristic  length  and  velocit:^  of  the  local  field,  both  of 
which  depend  upon  the  scale  of  the  region  used. 


'jsX-vi a ’ ■ Si 
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Appendix  A 


In  1.1  the  ten-  or  virial  equations  were  develooed  in 
cartesian  coordinates.  Lot  us  now  ver  briefly  outline  their 
aeveloonent  in  a general  coordinate  system.  Let  the  coordinates 
of  the  vth  particle  o j Let  the  ocsition  of  the  vth  particle 

relative  to  the  or^'p  . he  ^eoresemoo  oy  the  vector  q^)  . We 

note  that 


.j  _k> 


=§!iVJ 


The  Lagrengian  equation.,  of  r.hion  are 


. ■t/jffJ 1 _ = o1 

''bVV  t/  V 


and  are  readily  spm  to  be  covariant  and  of  rank  one  because  L 
is  a scalar,  dt  is  an  invariant,  and  tb  vQ^  satisfy  (1.1.2)  in 
which  8 VW  is  a scalar.  M 1 ti-ol^inp  by  (2)  may  be 

written  as 

£rgby)X|  . ?i . ♦ |j  (yk>  vo.  c 

L vq  J °vq 

(1.1.5)  — (1.1.8)  are  redefined  as 


J?i  - 


3 (,.qk)  -% 


3 = r sJ  Ak  _3L 


^■•k  ^ 3 51 

V i 


®fi  =rjbv,k) 


nfi  = n^j  <v^v<k 


Thus  (3)  is  the  mixed  tensor  equation  of  rank  two 


11  J.ji  = 2TiJ  + * ,3i  + 


wl  f’.V'  NV»?,.i|*r*ri  • <9  Arriba  4"  iorijB 

pl  = v1 

v v dt  v 


which  are  contr&variant , we  should  have  obtained  the  contravar iant 
form  of  (8),  in  which 


r*J  = 


vm^(vqk)  i?i(vqn) 


(10) 


2mij  = 


vm|i(vqk)  f:(vqn) 


<***  = ZZrlfiJ*) 


(12) 
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Appendix  B 


In  1.1  we  tacitly  assumed  that  the  dynamics  of  our  svstem  of 
particles  could  be  described  in  terms  of  the  snace  coordinates 


of  each  particle  and  the  time  der i,r; 


,,q' 


These  coord- 


inates do  not  include  spin  coordinates  for  the  practical  reason 
that  we  would  then  be  faced  with  working  in  a space  with  a much 


more  complicated  geometry.  It  follows,  then,  that  if  dipoles, 
quadripole?,  etc.  are  to  be  considered,  they  must  be  decomposed 
into  two,  four,  or  more  "elementary”  particles  with  suitable 
constraints . 


With  this  decomposition  it  follows  that  the  field  of  any 
particle  must  be  spherically  symmetric  in  at  least  the  proper 
coordinate  system  of  the  particle.  Hence  the  angular  momentum 
of  the  system  is  conserved  and  the  portion  of  due  to  the 

interaction  with  each  other  of  the  particles  in  the  region  is 
symmetric  in  its  indices.  We  have  denoted  the  symmetric  portion 
by  ^ in  (1.1.20).  Thus,  if  V ( r)  denotes  the 
potential  energy  of  the  interaction  of  the  vtn  and  T}th  particle, 
we  have 


We  define  to  be  the  potential  energy  of  the  interaction  of  the 
particles  with  themselves.  Thus 


(1) 


where 


T 
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V1  is  the  local  portion  of  V. 

If  we  define  as  the  local  portion  of  L,  i.e.  due  to  the 
interaction  of  the  particles  in  the  region  with  each  other,  then 
from  (1.1.7)  we  have  in  cartesian  coordinates  that 


.. 4 * ■ — 

j _ % 


:i  fh 

3 x 
v 


J 


Tlj  = -ZL 


3 x^ 


= -J> 


T],V 


xi 

vx  3_.  r 


T)V  I— 


V(_.r) 


~r  3 r 
tiv 


t}v  T)v 


8vxJ 


a^.v(««r)  (y*1  - ■nxl)vxi} 


ip  w 

3 r 
T)Vr 

(T,Vr) 

2 

- AT"  r 

( 

(VX1  " T, 

xi)(vxj  ”uxJ) 

2b  nv 

3 r ] 

ivr)2 

+ (vxl 

— X1)  ( X*5" 

ri  'v 

+ ^ j 

1 

(nvr) 

j 

1 

2 T)v' 

r\,v  ' 

. 3r,vTivr> 

( x1  — x1)  ( x^-  XJ") 

V T)  'v  T)  ‘ 

8„r 

?)V 

. rs 

i _ r ) 
r^v 

(2) 


We  see  that  as  required.  Let  us  diagonalize 

by  referring  it  to  its  principal  axes.  We  have,  then,  only  the 
dec  „ 

three  terms  f - Repeated  Greek  indices  do  not  imply  summation 
convention . Taa  is  the  contribution  to  the  potential t ensor  of 
the  forces  and  displacements  in  the  a direction. 


If  the 


TjVv(r)vr)  are  homogeneous  of  degree  n,  then 

1 2 


P^V-t 


s^7v  ,ivv(rivr) 


a a 
x — x 

v n_ 

TIV1* 


(3) 
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We  define  the  diagonal  tensor  of  the  second  rank  (‘7'°’)^  by 


Using  (1)  and  (4),  (3)  may  be  written 

nf°'a  = - n( Y a)2  Vx  (5) 

The  '7^~a  may  be  interpreted  in  a sense  as  the  direction 
cosines  of  the  eccentricity  of  the  region.  We  note  from  (5) 
that 

£Z^a)z  = i 

a-]. 

Dynamical  symmetry  cf  the  region  gives  = T"2  = tT*.  An 

elongation  of  the  region  in  the  {3  direction  makes  larger  than 

the  other  two  Y~ , etc* 


4? 
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